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ABSTRACT 

Low-frequency radio surveys are ideal for selecting orientation-independent sam- 
ples of extragalactic sources because the sample members are selected by virtue of 
their isotropic steep-spectrum extended emission. We use the new 7C Redshift Survey 
along with the brighter 3CRR and 6C samples to investigate the fraction of objects 
with observed broad emission lines - the 'quasar fraction' - as a function of redshift 
and of radio and narrow emission line luminosity. We find that the quasar fraction 
is more strongly dependent upon luminosity (both narrow line and radio) than it is 
on redshift. Above a narrow [Oil] emission line luminosity of log 10 (L[oirj/W) > 35 [or 

radio luminosity log 10 (Li5i/ W Hz -1 sr -1 ) > 26.5], the quasar fraction is virtually 
independent of redshift and luminosity; this is consistent with a simple unified scheme 
with an obscuring torus with a half-opening angle #trans ~ 53°. For objects with less 
luminous narrow lines, the quasar fraction is lower. We show that this is not due to 
the difficulty of detecting lower-luminosity broad emission lines in a less luminous, 
but otherwise similar, quasar population. We discuss evidence which supports at least 
two probable physical causes for the drop in quasar fraction at low luminosity: (i) a 
gradual decrease in #trans and/or a gradual increase in the fraction of lightly-reddened 
(0 < Ay ;$ 5) lines-of-sight with decreasing quasar luminosity; and (ii) the emergence 
of a distinct second population of low luminosity radio sources which, like M87, lack 
a well-fed quasar nucleus and may well lack a thick obscuring torus. 
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1 INTRODUCTION 

Orientation-based unified schemes for radio-loud quasars 
and powerful radio galaxies play a key role in our under- 
standing of these objects. Ever since the conception of the 
idea that powerful radio galaxies are simply quasars with 
their jet axes oriented away from our line of sight, such that 
the nuclear continuum and broad-line regions are obscured 
by a dusty torus or warped disc (Scheuer 1987; Barthel 
1989), evidence has been sought to prove that these uni- 
fication schemes are correct. Reviews of unification schemes 
for active galactic nuclei have been presented by Antonucci 
(1993) and Urry & Padovani (1995). Although the majority 
of observations are consistent with them, some observations 
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have been used to cast doubt on their viability over all ranges 
of radio luminosities and redshifts. 

The narrow line emission in radio sources is observed to 
be emitted largely from beyond the obscuring material (e.g. 
McCarthy, Spinrad & van Breugel 1995; Hes, Barthel & Fos- 
bury 1996) and therefore is independent of the jet axis orien- 
tation. Hence in this model one would expect radio galaxies 
and quasars to have similar narrow line luminosities. The 
strong positive correlation between the extended radio lumi- 
nosities and narrow line luminosities of radio sources (Baum 
& Heckman 1989; Rawlings et al. 1989; Willott et al. 1999) 
means that the samples of radio galaxies and quasars to be 
compared must be matched in extended radio luminosity. At 
low redshift (z < 0.8), there have been claims that quasars 
are observed to have [OIII] line luminosities a factor of 5- 
10 greater than radio galaxies of similar radio luminosities 
(Baum & Heckman 1989; Jackson & Browne 1990; Lawrence 
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1991), although the [Oil] luminosities ol radio galaxies and 
quasars at these redshifts are indistinguishable (Browne & 
Jackson 1992; Hes et al. 1996). These differences between 
[Oil] and [OIII] have been interpreted as being due to par- 
tial obscuration of [OIII] since its higher ionization potential 
means it is likely to be emitted from closer to the nucleus 
than [Oil]. However, Jackson & Rawlings (1997) have inves- 
tigated the [OIII] luminosities of z > 1 radio galaxies and 
quasars and find their distributions indistinguishable. Us- 
ing a combined 7C/3CRR dataset Willott et al. (1999) find 
that quasars have more luminous narrow lines than radio 
galaxies at intermediate radio luminosities [26 < log 10 (Li5i 
/ W Hz _1 sr _1 ) < 27], but they are similar at higher ra- 
dio luminosities; this result leads to different slopes for the 
narrow-line versus radio luminosity correlation for quasars 
and radio galaxies. 

Barthel (1989) showed that in the redshift range 0.5 < 
2 < 1.0 the linear size distributions and fraction of quasars 
in the 3CRR sample are consistent with all radio galaxies 
having their jet axes at an angle #trans > 44° from our line- 
of-sight and all the quasars having #trans < 44° , where #trans 
is the half-opening angle of the obscuring torus; a value of 
0trans = 45° has been adopted by many on the basis of this 
paper. However, at higher redshifts in the 3CRR sample, the 
fraction of quasars increases giving # t rans ~ 60° (e.g. Singal 
1993). Note that because of the tight luminosity-redshift 
correlation inherent in a single flux-limited sample, this ap- 
parent correlation with redshift may be due instead to a 
correlation with radio luminosity. Singal (1996) used sev- 
eral large samples with differing radio flux-density limits to 
find that the quasar fraction in radio samples declines with 
decreasing radio flux-densities. However, due to the corre- 
lation between the optical and radio luminosities of steep- 
spectrum quasars (Serjeant et al. 1998; Willott et al. 1998a), 
radio-fainter samples will contain optically-fainter quasars. 
Therefore the quasar fraction in faint samples might be un- 
derestimated if only quasars brighter than a certain optical 
magnitude limit are identified, which may well be the case 
in some of these fainter samples used by Singal. 

Using the new 7C Redshift Survey, a low-frequency ra- 
dio sample selected at a flux-density limit 25 x lower than 
the 3CRR sample with « 90% spectroscopic redshift com- 
pleteness, we have previously shown that the quasar fraction 
does not depend strongly upon radio luminosity or redshift 
for 1 < z < 3 (Willott et al. 1998b). Note that the low- 
frequency selection of the samples is crucial, because this 
ensures selection on extended radio flux which should be en- 
tirely isotropic. In this paper we use a larger sample to inves- 
tigate how the quasar fraction of radio sources depends upon 
redshift, radio luminosity and narrow emission line luminos- 
ity and we consider the implications for orientation-based 
unified schemes. In a companion paper (Blundell, Rawlings 
& Willott in prep.) we will investigate the constraints placed 
on unified schemes by the radio properties of quasars and 
radio galaxies from the combined 7C/6C/3CRR dataset. We 
assume throughout this paper that H Q = 50 km s _1 Mpc _1 
and q = 0.5. 



2 THE SAMPLE 

Our sample consists of three completely-identified samples 
selected at (similar) low radio frequencies. These samples 
have different flux limits to provide broad coverage of the ra- 
dio luminosity-redshift (Lisx —z) plane (see Figure [i]). This 
enables the separation of evolutionary- and luminosity- 
dependent effects. 

Our faintest sample is the 7C Redshift Survey which is 
briefly described here (see also Willott et al. 1998a; Willott 
et al. in prep.; Blundell et al. in prep; Lacy et al. 1999). The 
7C Redshift Survey includes all sources with flux-densities 
at 151 MHz 5*151 > 0.5 Jy in three selected regions of sky 
(7C-I, 7C-II and 7C-III with a total sky area of 0.022 sr). 
The sample contains 130 radio sources which have all been 
identified with an optical/near-IR counterpart. The spec- 
troscopic completeness is « 90% with most of the remaining 
sources having redshifts well-constrained by optical/near-IR 
photometry (Willott, Rawlings & Blundell, in prep.). 

The bright sample used is the 3CRR sample of Laing, 
Riley & Longair (1983), which has complete redshift infor- 
mation for all 173 sources, selected with Sirs > 10.9 Jy. 
3C 231 (M82) is excluded because it is a nearby starburst 
galaxy and not a radio-loud AGN. 3C 345 and 3C 454.3 are 
flat-spectrum quasars which are excluded on the grounds of 
Doppler-boosting of compact components being responsible 
for raising their fluxes above the selection limit. One quasar 
from the 7C sample (5C7.230) was excluded for the same 
reason. 

The intermediate flux-density sample used is a revision 
of the 6C sample of Eales (1985) (see Rawlings, Eales & Lacy 
in prep., hereafter REL, for details). The flux limits of this 
sample are 2.0 < S151 < 3.93 Jy and the sky area covered is 
0.103 sr. Only 2 of the 58 sources in this sample do not have 
redshifts determined from spectroscopy. One object is faint 
in the near-infrared (K > 19), so we take it as a galaxy at 
z = 2.0 in this paper. The other is relatively bright in the 
near-IR but very faint in the optical. Its red colour suggests 
a galaxy with a redshift in the range 0.8 < z < 2 and we 
assume z = 1.4 here. Full details of the sample, including 
optical spectra, will appear in REL. 

All the sources in two of the three 7C regions, 7C-I and 
7C-II, and in the 6C sample have been imaged in the near- 
infrared K-band. This has enabled the discovery of a few 
lightly reddened quasars, which would otherwise have been 
identified as galaxies on the basis of optical spectroscopy (see 
e.g. Willott et al. 1998a). In addition we have obtained near- 
infrared spectroscopy of some possible reddened quasars in 
these samples to search for broad emission lines (Willott et 
al. in prep.). On the basis of the available imaging and spec- 
troscopy, all sources have been classified into one of three 
categories: quasar, broad line radio galaxy (BLRG) or ra- 
dio galaxy. Given the relatively small number of reddened 
quasars found in the 6C, 7C-I and 7C-II samples, we do not 
believe the lack of near-infrared data for the 7C-III sample 
has caused a statistically-significant number of quasars to 
be mis-classified. Spectra of the quasars and BLRGs in 7C- 
I and 7C-II are presented with a detailed discussion of the 
classification scheme in Willott et al. (1998a). 

Our total sample comprises 356 sources (one object, 3C 
200, is common to both the 3CRR and 7C-II samples). We 
do not exclude compact symmetric objects (CSOs) from our 



© 0000 RAS, MNRAS 000, 000-000 



The quasar fraction in low-frequency selected complete samples 3 



29 



27 



O 



25 



• 

• ( 

c 

o • 

o US* 5 


• 

IP JD 

o 


1111,1111,1111 

> 

V 

□ o_ 

: O A 

CA«* O Ji A 

□ 

A □ 

^ □ _ 
▲ 

□ 

□ 


W ' 

CO* ■ 

<X.O' 


♦■ft " i 

&n. ■ 

A^ A * 


.□A 


g 


35(8) 
32(6) 


46(8) 
37(6) 


39(9) 


14(4) 





29 



28 



27 



26 



25 



1111,1111 

• 

• - 

G 

o o # 


1 

• 

IP JO 

' 


I 

□ o_ 

□ a 

OA«* * A 

A □ 

A 

□ 


g*o ■ 

~ "f ^ 

•Of ^ ^ 


U ' 


u <P A n 
♦ fa A 

o A 

- • 

□ 
A 

□ 


f 

ITT 1 1 I T T 1 


38 (9) ; 49 (8) 
44(8) ; 37(7) 


43(9) 


20(10); 

1 , 



0.0 0.5 1.0 1.5 2.0 2.5 3.0 3.5 0.0 0.5 1.0 1.5 2.0 2.5 3.0 3.5 
Redshift z Redshift z 



Figure 1. The radio luminosity-redshift plane for the 3CRR (circles), 6C (diamonds), 7CT and 7C-II (triangles) and 7C-III (squares) 
complete samples in which only the 323 objects deemed FRIIs are plotted (see Section 2). Filled symbols are quasars or broad line radio 
galaxies and open symbols are radio galaxies. We have binned the L151 — z plane in luminosity and redshift, of which 6 bins are well 
populated by sources from the complete samples. The percentage of quasars along with associated Poisson errors in each bin are shown 
in the bottom-right corner of the plot. The left plot shows all the FRII sources in the complete samples, whereas the right plot shows 
only those in the high emission-line luminosity category (see text for details). Note that there is one additional radio galaxy in the 7C-III 
sample which may be at z > 4 (Lacy et al. 1999). 



analysis as previous authors have (e.g. Barthel 1989, Singal 
1996), since we have no a priori reason to expect them not 
to be part of the unified schemes. In Section ^ we repeat 
our analysis excluding the CSOs to check this assumption. 
Sources with FRI radio structure are however excluded. The 
main reason for this was a practical one: most (24 out of 
33) of the FRIs are from the 3CRR sample, and the exist- 
ing spectrophotometry of these objects is inadequate for the 
quantitative investigations of this paper. However, as dis- 
cussed fully in Section 4.5, the lack of quasars with FRI ra- 
dio structure (e.g. Falcke, Gopal-Krishna & Biermann 1995), 
and the concentration of these objects at low L151 and low 
z (see Fig. |l|) makes it easy to investigate the effects of the 
exclusion of FRIs on our study of the quasar fraction. Al- 
though the radio structures of the quasars 5C6.264 and 3C 
48 are arguably FRI, we count them in this paper as FRIIs 
because of their high radio luminosities. Indeed all objects 
that could not unambiguously characterised as FRI or FRII 
(e.g. the 'DA' sources in the notation of Blundell et al. in 
prep.) are also counted as FRIIs in this study. Thus our total 
combined sample contains 323 objects deemed to be FRIIs. 



3 THE QUASAR FRACTION 

We now investigate the fraction of objects which show broad 
emission lines (referred to hereafter simply as the quasar 
fraction) in the complete radio samples. Broad line radio 
galaxies (broad line objects with Mb > —23) are counted 



as quasars in this analysis, because they are most likely 
to be weak quasars viewed within the torus opening angle 
(Laing et al. 1994; Hardcastle et al. 1998; Fig. 6 of Willott 
et al. 1998a). However, objects with scattered broad lines 
seen only in polarised light (4 cases) are counted as galax- 
ies, because the fact that the broad lines are not observed 
directly indicates that the nuclear regions are heavily ob- 
scured, most likely by the torus. Since spectropolarimetry 
is unavailable for most sources, there is clearly some danger 
that the split between weak and scattered-light quasars is 
imperfect. A further concern is the question of the classifi- 
cation of lightly-reddened broad-line objects in which it re- 
mains uncertain whether (a) the reddening is related to the 
torus or is caused by dust further from the nucleus, for exam- 
ple in the host galaxy of the quasar; and (b) whether, partic- 
ularly given the inhomogeneous spectrophotometric dataset 
for the 3CRR sample, all such objects have yet to be dis- 
covered. Taking J; Ay £5 as a working definition of 'light' 
reddening, we have chosen to count all such objects known 
to us as quasars. We will return to this important point in 
Section 4.2. 

In Fig. [l] (left), we plot low- frequency radio luminos- 
ity L151 against redshift z for the combined sample of FRII 
sources. The L151 — z plane has been binned and the per- 
centage of quasars in each bin shown, along with the associ- 
ated Poisson errors. The first thing to note is that there are 
very few quasars at log 10 (Lisi/ W Hz -1 sr _1 ) < 26.5. Ex- 
cluding this region, the quasar fraction appears to increase 
slightly as a function of both redshift and radio luminos- 
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ity, ranging from 0.3 to 0.5. However, with Poisson errors of 
a ~ 0.08 these differences can at best be called marginally- 
significant, and we will not consider them further. Note also 
that for the high-luminosity bins, the median luminosity in 
each bin changes with redshift, so the effects of redshift and 
luminosity have not been completely separated. 

Laing et al. (f994) proposed that by excluding low- 
excitation radio galaxies (LEGs; Hine & Longair f979) 
which have very weak or absent emission lines and low- 
ionization narrow lines, the quasar fraction in 3CRR was 
not a function of luminosity (or redshift). Their classifica- 
tion for LEGs was that they have [OI1I] to Ha ratios of 
< 0.2 and [Olll] equivalent widths of < 3 A. Since the radio 
galaxies in the combined samples here span a large range of 
redshift, Balmer and/or [Olll] lines are often not observed in 
optical spectra and a classification scheme such as this can- 
not be applied. Instead we separate sources into those with 
high and low [Oil] emission line luminosities. The division 
between these classes we adopt is log 10 (L[on]/W) = 35.1. 
Below this line luminosity there are very few broad line ob- 
jects in the complete samples (Willott et al. 1999). This cut 
in narrow emission line luminosity was chosen to be equiv- 
alent to a rest-frame [Oil] equivalent width of 10 A for a 

quasar with Mb = — 23^|. Note that this division puts more 
sources in the low- luminosity category than those typically 
classified as LEGs. 

On the right-hand side of Fig. [j] we plot only sources 
with [Oil] luminosities log 10 (L [O n]/W) > 35.1. The differ- 
ences between the quasar fractions in the bins is reduced 
somewhat here. Now we find that (with the exception of 
the poorly-populated low radio luminosity bin) the quasar 
fraction is 0.40 in all the bins at the la level. The largest 
change has occurred in the low-redshift, intermediate lumi- 
nosity bin which contains many weak-lined 3CRR galaxies. 
Note that the 1 < z < 2, intermediate luminosity bin has 
the lowest quasar fraction. A possible reason for this is that 
this bin contains most of the galaxies without lines in their 
optical spectra which may not all truly lie within this red- 
shift range. Considering all 216 high emission line luminosity 
sources plotted on the right-hand side of Fig. [j], we find a 
quasar fraction of 0.40 ± 0.03 implying a mean torus half- 
opening angle of 53° ±3° for the luminous population. There 
is likely to be quite a range of torus opening angles and the 
small error presented here on the mean angle should not be 
interpreted as the dispersion in opening angles present in 
the population. 

We have repeated the analysis of this section excluding 
CSO sources (projected linear sizes < 30 kpc). We find that 

t The mean rest-frame equivalent width of the [OIII] line for BQS 
quasars in the study of Miller et al. (1992) is 30 A. Assuming a 
ratio of 3:1 for the [OIII]:[OII] equivalent widths, this translates 
to 10 A for [Oil]. This value was used in Willott et al. (1999) to 
determine the strength of the photoionizing continuum in radio 
sources. However, the typical ratio between [Oil] and [OIII] may 
actually be smaller than this, causing Willott et al to have over- 
estimated the mean [Oil] equivalent width by a factor of a few. 
Baker et al. (1999) finds a median [Oil] equivalent width of 4 A in 
the Molongo Quasar Sample. Also, in the LBQS composite quasar 
spectrum of Francis et al. (1991) the [Oil] equivalent width is 2 
A, which would give a line luminosity of log 10 -L[ou] = 35.1 for a 
Mb = —25 quasar. 



this reduces the quasar fraction by ~ 0.04 in all the bins. 
The quasar fraction of all luminous CSO sources is 0.56 ± 
0.08, different at the 2a level from that of non-CSO sources. 
We defer discussion of possible reasons for a higher quasar 
fraction for CSOs in low-frequency selected samples to a 
future paper. 

4 POSSIBLE CAUSES OF THE CHANGE IN 
QUASAR FRACTION WITH LUMINOSITY 

4.1 Selection effects caused by lower intrinsic 
broad line luminosity 

We first investigate whether the decrease in quasar fraction 
with decreasing narrow-line and radio luminosity is a sim- 
ple selection effect. If the low-luminosity objects have lower 
luminosity narrow emission lines it is natural to expect that 
their broad line and continuum luminosities should be lower 
too (e.g. Miller et al. 1992). It follows that their spectra 
may on average have lower signal-to-noise and weak broad 
lines may be missed. Given the strong correlation between 
narrow-line and radio luminosities (e.g. Willott et al. 1999) 
this could explain the deficit of broad line objects at radio 
luminosities log 10 (Li5i / WHz~ 1 sr~ 1 ) < 26.5, and also the 
apparent constancy of the quasar fraction for objects with 
luminous narrow lines discussed in Section 3. 

To test whether selection effects such as these could 
cause the apparent lack of low-luminosity broad line ob- 
jects, we first calculate the expected broad line fluxes of 
objects in the 3CRR and 7C samples as a function of lumi- 
nosity. To do this we make use of the fact that the narrow 
emission line luminosities of radio galaxies and quasars are 
positively correlated with low-frequency radio luminosity as 
shown in Willott et al. (1999). This correlation has a slope of 
0.79 ±0.04 in the sense that Lroii] oc L°'Ji and is most likely 
due to a correlation between jet power and photoionizing 
continuum luminosity (Rawlings & Saunders 1991). Use of 
this relation enables one to estimate the narrow line [OR] 
luminosity for values of radio luminosity. 

A similar correlation between the broad line luminosity 
and radio luminosity is expected to hold. There are several 
lines of evidence supporting this. First, radio luminosity is 
correlated with the optical continuum luminosity for steep- 
spectrum quasars (Serjeant et al. 1998) and quasars show 
little variation in broad line equivalent widths over a wide 
range of luminosity (e.g. Osmer & Shields 1999). Second, 
Celotti et al. (1997) have shown that broad line luminosities 
are correlated with the power in pc-scale jets in radio-loud 
quasars. 

To determine the relationship between the strengths of 
narrow and broad lines in quasar/BLRG spectra, we adopt 
the line ratios from the composite quasar spectrum of Fran- 
cis et al. (1991). Although this is from an optically-selected 
sample, we do not expect significant differences between 
the line ratios of radio-loud and radio-quiet quasars. This 
is borne out by a comparison with the Molonglo Quasar 
Sample of Baker & Hunstead (1995): repeating the analysis 
using their line ratios gives a virtually identical result. 

We consider a range of radio luminosities and attach to 
each one typical redshifts at which sources of these luminosi- 
ties are observed in both the 3CRR and 7C samples. The cor- 
relation between the [Oil] line and radio luminosity is then 
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Table 1. This table shows typical redshifts of sources from the 
3CRR and 7C samples (columns 2 and 3) at various radio lumi- 
nosities (column 1). Columns 4 and 5 show the strongest broad 
lines observable over the relevant optical wavelength ranges cor- 
responding to these redshifts. 

used to determine typical [Oil] luminosities at each value of 
radio luminosity. From the typical wavelength range covered 
by optical spectra (4000 A- 8500A), we find the strongest 
broad line observable at each redshift. The typical luminosi- 
ties, redshifts and strongest broad lines are shown in Table 
[j] Finally, we use the ratio of broad to [Oil] lines in Francis 
et al. (1991) to calculate the luminosity and hence flux of 
the strongest broad line expected at each luminosity-redshift 
pair. 

There is considerable scatter in the relationship be- 
tween narrow emission line luminosity and radio luminosity 
(a = 0.5 dex). Extra scatter comes in due to the conversion 
from narrow to broad line luminosities. Therefore a total 
scatter of la = 0.6 dex was assumed to determine the frac- 
tion of objects with broad line fluxes significantly below the 
characteristic values. 

Figure ^ shows the expected flux of the strongest broad 
line observable at each redshift as listed in Table hi for both 
3CRR and 7C. Also plotted are the actual data for the 7C-I 
and 7C-II quasars and BLRGs (Willott et al. 1998a) and all 
3CRR quasars with z < 0.86 in the RA range 18-12 hr (all 
from Jackson & Browne 1991 except 3C 109 from Goodrich 
& Cohen 1992). In general points corresponding to 3CRR 
BLRGs are not plotted because only a few of their broad line 
fluxes are available in the literature, although a few points 
from Hill et al. (1996) are plotted to illustrate the discussion 
of Section 4.2. The quasars which are red(dened) [a op t > 1, 
where /„ oc v~° L ° vt ] are shown as open symbols on Fig. ^; 
again the 3CRR spectrophotometric dataset is too inhomo- 
geneous for this to be done reliably for 3CRR sources. Note 
that very few quasars fall below the limits of the model error 
bars. Virtually all of those which do are reddened (the ex- 
ception is 5C6.282 at z ~ 2, which has an unusual spectrum 
with very narrow Lya; see Willott et al. 1998a). The sim- 
ilarity of the predicted fluxes with those observed justifies 
the use of the line ratios and method described here. 

The key point to note from Fig. |^ is that as one consid- 
ers lower radio luminosities in each sample (i.e. lower red- 
shifts), the expected broad line fluxes increase. Therefore, 
so long as the spectra of lower redshift sources have similar 
exposure times to those at higher redshift, the broad lines 
should be visible if they are there and unreddened. For the 
3CRR sample the spectra are rather inhomogeneous and it 
is difficult to prove that this is the case (but note that the 




Figure 2. Expected broad line fluxes for 'typical' 3CRR (large 
circles) and 7C (large triangles) sources as a function of redshift 
according to the radio-optical correlation as described in Section 
4.1. The error bars show approximate la deviations (i.e. ?s 1 in 
6 sources should fall below the bottom of the error bars). The 
small symbols show the actual data for 3CRR (circles) and 7C 
(triangles) quasars and the 7C BLRGs (asterisks). Open triangles 
are 7C quasars which are known to be reddened. The open circles 
show the reddened quasar 3C 22 at z = 0.93 (broad Mgll upper 
limit from Rawlings et al. 1995) and the 3C BLRGs from Hill et 
al. 1996 (see Section 4.2). 

deficit of 3CRR sources is most marked at z < 0.3, where 
the Ha line is covered by optical spectra and the expected 
broad Ha fluxes are high) . However, in the 7C sample, simi- 
lar exposure times have been attained for the spectra at low 
redshift to those at high redshift. Therefore we are confident 
that in the absence of significant reddening of the BLR, an 
insignificant number of sources would be mis-classified. 

A further problem is whether a weak quasar spectrum 
can be discerned against a stronger host galaxy spectrum. 
To determine if this could cause broad line selection prob- 
lems, synthetic spectra of quasars and galaxies were created 
and combined. The model used for the galaxy spectra was 
a 1 Gyr old stellar population synthesis model of Bruzual & 
Chariot (1993). For the quasar spectra, the LBQS composite 
of Francis et al. (1991) was used. Poisson noise was added 
and it was found that broad Ha or Mgll should still be 
clearly visible for a quasar 2 magnitudes (in B-band) fainter 
than the galaxy (Fig. ^). Broad H/3, however, is more diffi- 
cult to detect. Spectra with poor blue wavelength coverage 
(e.g. no data below 5000 A) of sources at 0.3 < z < 0.8 would 
not include Ha or Mgll, in which case H/3 is the brightest 
line. Almost all the 7C and 6C spectra cover a large enough 
wavelength range to avoid this problem. However many spec- 
tra of 3CRR sources in the literature have a smaller wave- 
length range. But the radio-optical correlation (and Table ^) 
show that quasars at these redshifts in 3CRR typically have 
high luminosities and should not be fainter (at B-band) than 
their host galaxies. Indeed, the deficit of broad line objects 
in 3CRR only really begins at z ^ 0.3. Hence quasar broad 
lines should be clearly visible against typical background 
galaxy spectra down to log 10 (Li5i / WHz _1 sr _1 ) w 25.5, 
with possible problems only at fainter radio luminosities. 

Another possibility to consider is that the narrow lines 
may be much stronger than the broad lines, so only weak 
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Figure 3. Superposition of a blue model quasar/BLRG spectrum 
and a red galaxy spectrum as described in the text. The quasar 
is two magnitudes fainter in B— band than the galaxy. The upper 
line shows the sum of the two components. Note that the broad 
Mgll and Ho lines are clearly visible in the total flux spectrum, 
the former because of the quasar dominance at wavelengths A < 
3000 A and the latter because of its high equivalent width. 

wings in the line profiles would be seen. However, this pic- 
ture seems unlikely for many sources in the absence of broad 
line reddening, since if both the narrow and broad lines are 
photoionizied by the same nuclear source it would require a 
conspiracy between broad and narrow line covering factors. 

We conclude that the low quasar fraction at low lumi- 
nosity is not due to selection effects caused by the difficulty 
of detecting intrinsically less luminous broad lines. 

4.2 Increased reddening at low luminosities 

In this section we consider the possibility that the drop in 
quasar fraction at low luminosities is due, at least in part, to 
an increase in the fraction of lightly-reddened (0 ^ Ay < 5) 
lines-of-sight as quasar luminosity decreases. It can be seen 
from Fig. m that the majority of the 7C quasars are lightly 
reddened at the highest redshifts. We exclude from this pa- 
per any discussion of any systematic changes of quasar red- 
dening with cosmic epoch. 

The question of the range of reddening towards the 
nuclei of powerful (3CRR) z ~ 1 radio sources has re- 
cently been addressed by the 3.5/im imaging programmes of 
Simpson, Rawlings & Lacy (1999) and Simpson & Rawlings 
(2000). They concluded that these high luminosity nuclei are 
mostly either naked (Av ~ 0) or heavily obscured (Av ^ 15) 
with only a small fraction (~ 15 per cent) of intermediate 
lightly-reddened cases. If these results could be extrapolated 
to the lower narrow-line and radio luminosity regimes then 
light reddening of quasar nuclei would not be a serious issue 
for our study of the quasar fraction. 

However, two studies suggest that the fraction of 
lightly-reddened lines-of-sight might increase significantly 
at lower narrow- line and radio luminosities. First, Hill, 
Goodrich & DePoy (1996) studied the Paa and Ha lines 
in a complete sample of low redshift (0.1 < z < 0.2) 3CR 
sources, looking for evidence of reddened broad lines in ra- 
dio galaxies. These objects lie near the top of lowest-Lisi 



(and z) bin in Fig. m They found broad lines in six ob- 
jects out of a complete sample of thirteen (including two 
FRI sources), of which three had intermediate values of Ay- 
This study, although subject to a major problem with small 
number statistics together with the difficulties of detecting 
lightly reddened broad lines over certain ranges of redshift 
(see Fig. ^), hints that an increase in the fraction of the 
lightly-reddened population at low luminosities might pro- 
vide part of the explanation for the dropping quasar frac- 
tion at low luminosities. Second, Baker (1997) showed that a 
typical low-luminosity lobe-dominated quasar, i.e. an object 
akin to those near the top of lowest-Lisi (and z) bin in Fig. 
[l], has a reddening Av ~ 3 and is thus lightly reddened. 

We conclude that some of the drop in the quasar frac- 
tion at low luminosities might be explained by an increased 
chance of light reddening towards low luminosity quasar nu- 
clei. In other words, we cannot allay the suspicion that some 
lightly reddened quasars have been misclassified as galaxies 
on the basis of existing, chiefly optical, spectroscopy. In the 
lowest-Lisi and z bin of Fig. [I] this is a problem for (very 
low-z) 3CRR sources because of the inhomogenous optical 
dataset; it is a problem for the 6C and 7C sources because 
their redshifts are sufficient to remove the strong Ha line 
from the optical window (see Fig. ^ . To pursue this idea fur- 
ther we next consider specific physical models which might 
explain this behaviour. 



4.3 The receding torus model 

Lawrence (1991) proposed that the apparent increase in 
quasar fraction and hence # tra ns with luminosity in the 
3CRR sample could be explained naturally by a 'receding 
torus' model. Dust sublimates at T ~ 1500 K and there- 
fore the inner radius of a dusty torus/disc may be governed 
by the radius at which this temperature is achieved via ra- 
diation from the AGN. The more luminous the AGN the 
larger this radius and if the scale height of the torus h is 
independent of luminosity, then more luminous sources will 
have a larger 9 tI a,ns; the inner radius of the torus r scales 
as L 0,5 and Otmns = tan -1 (r/h). This model can also ex- 
plain more lightly-reddened quasars at low luminosity since 
a larger fraction of the lines-of-sight pass through interme- 
diate columns of obscuring material when the inner wall of 
the torus lies closer to the nucleus (see e.g. Fig. 8 of Hill et 
al. 1996). 

Simpson (1998) has shown that if the height of the torus 
is independent of the ionizing luminosity L then the receding 
torus model predicts a theoretical quasar fraction / q given 
by 

/q = l- ( 1+ T^ tan2 ^) 05 ' W 

where Lo and 6o are the normalisation luminosity and the 
torus half-opening angle at that luminosity, both of which 
are fixed by the measurement of the quasar fraction at some 
normalisation luminosity. We further assume that the [Oil] 
luminosity scales with the ionizing continuum luminosity, 
and that it is not significantly affected by the changing torus 
opening angle, so that we can replace ionising luminosities 
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Figure 4. The points with Poisson error bars show the quasar 
fractions in four bins of [Oil] emission line luminosity plotted 
against the median log 10 (L[on]) °f the sources in each bin. The 
data used were for the objects deemed FRIIs in the combined 
sample 3CRR, 6C and 7C dataset (see Willott et al. 1999). The 
line is a simple prediction made on the basis of the receding torus 
model, normalised at logio(L[oii]/W) = 35.5, as described in 
Section 4.3. 

in Eqn ^ by values of Ljaii]. ^ Binning all FRII sources in 
terms of I/[on] we find that for 35 < log 10 (L[on]/W) < 36, 
the quasar fraction is 0.32±0.04 giving 9q = 46° . This is used 
to normalise Equation |l| which is plotted in Fig. ^ along with 
the quasar fraction measured in four bins of log 10 (L[on])- 

This simple model appears to fit the data fairly well 
although it does predict a marginally sharper transition be- 
tween low and high quasar fractions than the data (at about 
the 1<t level). Note that the above model depends crucially 
upon the height of the torus being independent of luminos- 
ity. If the height of the torus increases with luminosity, albeit 
more slowly than the inner radius, then this would flatten 
the model distribution to something more closely resembling 
the data. 

Taken together with the results of Hill et al. (1996), we 
conclude that the receding torus model provides a plausible 
if non-unique explanation for the declining quasar fraction 
below log 10 (L [O n]/W) = 35.1. 

4.4 Dust— ejection by jets or outflow 

Tadhunter et al. (1999) have recently obtained high- 
resolution near-infrared imaging of the powerful radio galaxy 
Cygnus A which shows an edge-brightened biconical struc- 
ture roughly aligned with the jets. They suggest that at 

■f As emphasised by Simpson, [Oil] is probably a far from ideal 
choice of narrow emission line for this procedure because for an 
ionisation-boundcd emission line cloud, the [Oil] emission is from 
a region whose size depends much less strongly on the incident 
photoionising flux than the sizes of the regions responsible for 
higher-ionization lines like [OIII]. This leads to a much weaker de- 
pendence of [Oil] luminosity on incident photoionising luminosity 
L than the near proportionality between, say, [OIII] luminosity 
and L. Unfortunately, the choice of [Oil] was forced on us by the 
available spectrophotometric data (see Willott et al. 1999). 



least some of the very high reddening towards the nucleus 
of Cygnus A (Av ~ 70 mag) occurs in a kiloparsec scale 
dust-lane. For radio-loud sources, the jets may clear a path 
through this dust, so that only when viewed within a crit- 
ical angle of the jet-axis is there a chance of a line-of-sight 
having little reddening. The jet power could control the dust- 
clcaring efficiency in such a way as to give an anti-correlation 
between reddening and radio luminosity. The conceptual dif- 
ference between this model and the receding torus model is 
that it is now a decrease in jet power (and hence the radio 
luminosity) rather than a decrease in optical luminosity that 
controls the closing down of the effective opening angle of the 
radiation cone at low luminosities. Tadhunter et al. specu- 
late that lightly reddened quasars are objects viewed within 
the nuclear torus opening angle, but not within the region 
cleared of dust by the jets. Another possibility (e.g. Dopita 
et al. 1998) is that dust is cleared by a radiation-driven out- 
flow in which the opening angle and strength of the outflow 
could be related to the luminosity of the nucleus and hence 
give rise to the observed dependence of quasar fraction on 
luminosity. 

4.5 Dual-population models 

The evidence of increased reddening at low luminosities 
(Section 4.2) and the quantitative success of the receding 
torus model (Section 4.3) hint that we have a reasonable ex- 
planation for at least part of the drop in quasar fraction at 
low luminosities. In this section we discuss separate evidence 
for a second conceptually different mechanism. 

A corollary of Fig. ^ is that at low luminosities the low 
fraction of quasars results from a small value for the typical 
opening angle of the torus #trans- Falcke et al. (1995) have 
cited this as a strength of the receding torus model since if 
the opening angle approaches the beaming angle of the jet 
then there will be an unobscured view of the nucleus of low 
power jets only when their emission is also strongly Doppler- 
boosted. Such objects would be classified as BL Lac objects 
rather than quasars, providing a neat explanation for the 
well-known lack of quasars with FRI radio structures. We 
are unaware of any systematic measurements of broad line 
fluxes in BL Lac type objects which would provide a quanti- 
tative test of this hypothesis. However, recent observations 
of BL Lacertae itself (Corbett et al. 2000), an object with 
an extended radio flux characteristic of a low-luminosity ra- 
dio galaxy (Antonucci & Ulvestad 1985), suggest that in at 
least this object there are broad emission lines at about the 
predicted level, and that these lines are excited by a hot 
accretion disc. 

However, other recent results suggest that the small 
opening angles required at low luminosities from Fig. ^| are 
more likely to argue against the generality of the receding 
torus model. If it were true that low luminosity radio galax- 
ies possessed very thick obscuring tori then it should prove 
very difficult to get a clear view of their central regions at op- 
tical wavelengths. Chiaberge, Capetti & Celotti (1999) have 
studied the nuclear emission from FRI galaxies using the 
HST. In most cases they find probable optical synchrotron 
from the nuclear regions, which they interpret as implying 
a clear view of the inner jet and hence a lack of obscuring 
material on scales greater than ~ 100 pc. This straight- 
forwardly rules out a kpc-scale obscuration model for FRIs 
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along the lines of the model proposed by Tadhunter et al. 
(1999) for Cyg A. Still stronger constraints come from the 
nearby FRI radio galaxy M87 since its optical synchrotron 
core is smaller than ~ 5 pc and varies on the timescales 
of months (Tsvetanov, Kriss & Ford 1997), implying that 
the emission region is a pc or smaller in size. Since the in- 
ner regions of the optical synchrotron-emitting jet are in 
clear view, and these size scales are comparable in size to 
the expected size of the obscuring torus in low-luminosity 
quasars (e.g. Netzer & Laor 1993), these results begin to 
question the existence of any thick torus, at least in the 
case of M87. Chiaberge et al. maker stronger statements 
than this by making the additional assumption that this op- 
tical emission is emitted from within the region probed by 
radio VLBI which is argued to be within ~ 0.1 pc of the cen- 
tral black hole (Junor, Biretta & Livio 1999). We remain to 
be convinced that these observations rule out the existence 
of a pc-scaled torus in M87. There is a plethora of evidence 
from studies of M87 using optical continuum (Chiaberge et 
al. 1999), emission-line (Ford et al. 1994) and hard X-ray 
(Guainazzi & Molendi 1999) observations that it contains a 
supermassive black hole accreting at a tiny (< 10 -4 ) frac- 
tion of its Eddington-limited rate (see Willott et al. 1999 for 
a discussion of this object in the context of the radio-optical 
correlation) . 

As shown in Fig. [I], excluding objects with weak nar- 
row emission lines gives a quasar fraction which is virtu- 
ally independent of redshift and/or radio luminosity. This 
observation can be naturally explained if the weak-line ob- 
jects belong to a different population from the strong emis- 
sion line objects. Perhaps, as is seemingly the case for M87, 
a large fraction of these objects are fundamentally differ- 
ent in that they lack a well-fed quasar nucleus and also a 
thick obscuring torus (Chiaberge et al. 1999), and perhaps 
their jets are powered by a fundamentally different mech- 
anism (e.g. Blandford & Begelman 1999). The continuity 
of the radio-optical correlation for radio sources provides a 
weak argument against a fundamental change in the accre- 
tion process (Rawlings & Saunders 1991), but evidence for 
systematic changes in the slope of this relation at low lumi- 
nosities (Zirbel & Baum 1995) may be evidence that there 
are at least some differences at low radio luminosities. Argu- 
ing along similar lines, Laing et al. (1994) proposed that the 
low-excitation radio galaxies (LEGs) in their 3CRR sample 
would lack broad lines if viewed from any orientation. Re- 
call that we specifically excluded FRI radio sources from our 
analysis of the quasar fraction in Section 3. The only differ- 
ence that including the FRIs would have made, would be 
to decrease the quasar fraction in the lowest luminosity bin 
from 0.14 ± 0.04 to 0.10 ± 0.03. 

Such a 'two population' model would remove another 
troubling problem for any extrapolation of the receding 
torus model to the lowest quasar luminosities. Seyfert galax- 
ies - radio-quiet AGN with similar optical narrow-line lu- 
minosities to radio galaxies in the lowest-Li5i bin of Fig. 

- do not have the narrow opening angles predicted by 
the receding torus model. Direct observations of the ioniza- 
tion cones in nearby Seyfert galaxies typically show cones 
with half-opening angles in the range 20 — 50° (Wilson & 
Tsvetanov 1994), considerably greater than those predicted. 
Moreover, the vast majority of Seyfert galaxies do not seem 
to be obscured, i.e. of type 2, as would be required by nar- 



row opening angles. Hence, as noted by Falcke et al (1995) 
it would appear that the receding torus model cannot ap- 
ply similarly to radio-quiet and radio-loud AGN requiring 
some ad-hoc explanation, for example a different geometry 
for the obscuring material (i.e. a smaller scale height), due to 
a difference in black hole mass, environment and/or angular 
momentum. 

We can make a crude association of the putative sec- 
ond population of radio sources with all sources with nar- 
row emission line luminosities log 10 (L[on]/W) < 35.1. Due 
to the scatter in the radio-optical correlation, it is not ex- 
pected that this would lead to a clean division between the 
two populations at a particular [Oil] luminosity, and if the 
data were available, a classification based on line ratios, and 
hence excitation, might well prove cleaner. However, sub- 
ject to this limitation, the change in quasar fraction with 
redshift of a combination of the two populations (seen on 
the left plot of Fig. |l|) is then naturally explained by the 
less rapid cosmic evolution of the low-luminosity population 
(e.g. Urry & Padovani 1995) than the high-luminosity pop- 
ulation (shown on the right plot of Fig. hi). Hence at high- 
redshift, the contribution of the low-luminosity population 
is negligible and the quasar fraction of 0.4 is just that of 
the high-luminosity population. Dual population modelling 
of the low-frequency radio luminosity function is consistent 
with just such a scheme (Jackson & Wall 1999; Willott et 
al. in prep.). Also, Laing et al. (1994) and Hardcastle et 
al. (1998) find that low- excitation radio galaxies have linear 
size and core prominence distributions consistent with an 
isotropic population. 

There are two residual concerns with the two population 
model. First, many of the low-luminosity objects have emis- 
sion lines, so their excitation and the correlation between 
luminosity and ionization parameter (Saunders et al. 1989; 
Tadhunter et al. 1998) need some explanation. Second, it 
must explain why radio galaxies and quasars have different 
narrow line luminosity distributions at intermediate lumi- 
nosities (Jackson & Browne 1990) but similar distributions 
at high luminosities (Jackson & Rawlings 1997). The first 
concern can be addressed by simple analogy with M87: the 
absence of a broad-line quasar nucleus does not mean an 
absence of a photoionising source; for example Dopita et al. 
(1997) favour radiative shocks as a source of the excitation 
for the Ha lines in the nuclear disc of M87, and shock models 
in which ionization parameter correlates with line luminos- 
ity are easily envisaged (e.g. Dopita & Sutherland 1995). 
The second concern is probably also easily dealt with. Con- 
sidering first the 3CRR sample, the dual-population model 
is consistent with the drop in quasar fraction at low lumi- 
nosities because a low luminosity 3CRR source is necessar- 
ily at low redshift where there is clearly a mixture of both 
populations; at high redshifts only the high-luminosity pop- 
ulation is observed. Thus any comparative narrow emission 
line study of quasars and radio galaxies which is based on 
bright radio samples (Jackson & Browne 1990; Jackson & 
Rawlings 1997) should yield different line luminosity dis- 
tributions at low redshift, and similar distributions at high 
redshift, which is just as observed. However, a small dif- 
ference in the distributions of the [Oil] line luminosities of 
intermediate luminosity quasars and radio galaxies is also 
observed in the 7C sample (Willott et al. 1999). The lower 
radio flux limit of this sample means that these intermediate 
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luminosity sources are at high redshift, where there are few 
low emission line luminosity objects, so the difference is not 
due to the mixing of the two populations. This is most likely 
the result of small but inevitable biases pointed out by Rawl- 
ings & Saunders (1991), and quantified in the context of the 
receding torus model by Simpson (1998): a positive correla- 
tion between quasar luminosity and opening angle, coupled 
with inevitable scatter in quasar luminosity at a fixed radio 
luminosity, means that objects viewed within the opening 
angle are biased towards the more luminous objects within 
the scatter. 

4.6 Time variability 

There is yet one more possible cause of the drop in quasar 
fraction at low luminosities: systematic differences in the 
time variability of objects with luminosity. 

In Willott et al. (1999) we used the narrow emission 
line-radio correlation to suggest that the most luminous ob- 
jects are probably accreting at rates close to the Edding- 
ton limit, but lower luminosity sources are sub-Eddington 
accreters. It therefore seems plausible that the lower lumi- 
nosity sources have more scope for variability, since an ob- 
ject accreting at the Eddington rate should have a ready 
fuel supply which is accreted at a fairly steady rate lim- 
ited by radiation pressure. In contrast, sub-Eddington ac- 
cretion suggests there is not a ready supply of fuel available 
and hence fluctuations in accretion rate may be more likely. 
Indeed observations of radio-quiet quasars show that the 
lower luminosity quasars are more highly optically-variable 
than higher luminosity quasars over timescales of a few years 
(e.g. Veron & Hawkins 1995; Cristiani et al. 1996; Paltani 
& Courvoisier 1997) , although it should be noted that some 
authors attribute at least some of this variability to gravi- 
tational microlensing (e.g. Hawkins & Taylor 1997). 

The small quasar fraction at low luminosities could be 
explained by these objects spending a significant fraction 
of their active lifetimes in a 'quiet' state whereas the more 
luminous objects are continuously active over their entire 
lifetime. Due to light travel time effects, different emission 
regions of quasars have different variability timescales: BLR 
~ months; NLR ~ 10 4 yr; radio lobes ~ 10 6 yr. Reverbera- 
tion mapping of quasars has shown that the broad line fluxes 
follow the nuclear continuum flux with just such a time lag 
(see Peterson 1993 for a review) . If a quasar undergoes high- 
amplitude variability over a timescale ~ 100 yr, then only 
the continuum and BLR fluxes would be observed to un- 
dergo this strong variability, and the NLR and extended 
radio emission would simply reflect the time-averaged out- 
put of the central engine. Note that a quasar undergoing 
a decrease in luminosity of this sort of timescale may then 
appear as a low-excitation radio galaxy. If this behaviour 
occurs only in a significant fraction of the intermediate and 
low luminosity objects, then it could reproduce the observed 
change in quasar fraction. 

Essentially this model is a slightly more complicated 
version of the two-population model considered in Section 
4.5, with the two populations being unified along a time 
rather than a jet orientation axis. It is important to stress 
that there is as yet no firm evidence for the type of time 
variability required. However, the few long-term time vari- 
ability studies hint that there may be some interesting ef- 



fects: the BLRG 3C 390.3 was found to undergo a sustained 
decrease of 1.5 mags in the optical over « 80 years (Can- 
non, Penston & Penston 1968); however, Angione (1973) 
used archival plates to study the variability of 23 quasars 
over w 60 years finding no evidence for long-term sustained 
increases or decreases. 

A separate issue related to time variability is whether 
the probability of individual radio sources being observed to 
be quasars varies systematically throughout their lifetimes. 
Because flux-limited samples introduce inevitable biases into 
the age distributions of the sources they contain (e.g. Blun- 
dell, Rawlings & Willott 1999), this can lead to subtle ef- 
fects. 



5 CONCLUSIONS 

Using complete, low-frequency selected samples of radio- 
loud AGN we have investigated the fraction of observed 
broad line objects - the quasar fraction - as a function of 
redshift, and radio and emission line luminosity. Our find- 
ings are interpreted in terms of orientation-based unified 
schemes. We find that 

• Considering only those sources with strong narrow 
emission lines [log 10 (L[on]/W) > 35.1, corresponding 
roughly to log 10 (Li5i/W) > 26.5], the quasar fraction is 
virtually independent of radio luminosity and redshift, and 
is consistent with a simple unified scheme involving an ob- 
scuring torus with a half-opening angle 8 tra ,ns of 53° . 

• For less luminous emission-line sources, the quasar frac- 
tion is much lower, a finding which is not simply the result 
of selection effects induced by the lower intrinsic broad line 
luminosities of these sources. 

We have found evidence which supports at least two 
probable physical causes for the drop in quasar fraction at 
low luminosity: (i) a gradual decrease in # t rans and/or a grad- 
ual increase in the fraction of lightly-reddened (0 ^ Ay ^ 5) 
lines-of-sight with decreasing quasar luminosity; and (ii) the 
emergence of a second population of low luminosity radio 
sources which, like M87, lack a well-fed quasar nucleus, and 
may well lack a thick obscuring torus. 
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